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Abstract 

Nonlinear optical signals from an assembly of N noninteracting particles consist of an incoherent 
and a coherent component, whose magnitudes scale ~ N and ~ N{N — 1), respectively. A unified 
microscopic description of both types of signals is developed using a quantum electrodynamical 
(QED) treatment of the optical fields. Closed nonequilibrium Green's function expressions are 
derived that incorporate both stimulated and spontaneous processes. General (n + l)-wave mixing 
experiments are discussed as an example of spontaneously generated signals. When performed 
on a single particle, such signals cannot be expressed in terms of the nth order polarization, as 
predicted by the semiclassical theory. Stimulated processes are shown to be purely incoherent in 
nature. Within the QED framework, heterodyne-detected wave mixing signals are simply viewed as 
incoherent stimulated emission, whereas homodyne signals are generated by coherent spontaneous 
emission. 

PACS numbers: 42.65-k, 12.20.-m, 42.62. Fi, 42.50.-p, 42.50.Ct 
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I. INTRODUCTION 



Nonlinear optical processes in bulk materials are traditionally classified as either coherent 
or incoherent l]-^]. This classification refers to the relation of the macroscopic signal 
intensity of the sample to the nonlinear response of the individual molecules. Processes 
where the detected signal is obtained by simply adding up the contributions of the individual 
particles in the sample, are termed incoherent. Such signals hence scale as the number of 
particles N in the interaction volume. With coherent processes, on the other hand, the 
contributions of the constituent particles add up on the level of amplitudes rather than 
intensities, and the detected signal is proportional to iV 2 . 

Coherent optical signals are usually calculated semiclassically by adopting a two step 
procedure 3] - {4]. First a (linear or nonlinear) polarization P^ induced in the system 
by interactions with classical external fields is calculated microscopically. In the second 
step, P^ is used as a source in the macroscopic Maxwell equations to generate the signal. 
Frequency-domain signals are given by the absolute square of an amplitude related to the 
susceptibilities which, in turn, are proportional to N. This gives rise to the ~ N 2 
dependence of coherent signals Examples are Rayleigh scattering (n = 1), sum 

frequency generation (n = 2), four-wave mixing (n = 3). In contrast, signals such as Raman, 
fluorescence, two photon fluorescence etc. are incoherent and may not be calculated by the 
above semiclassical approach. Instead, they require a quantum description of the field to 
account for spontaneous emission. 

Treating both types of processes by a unified approach is of fundamental interest [lj], 
particularly when they coexist. For example, two photon fluorescence and second harmonic 
generation have been observed in the same sample p|. There is a considerable current 
effort to perform nonlinear optical measurements on small samples and even single nano 
particle, or m „>ecu>es fl-fl Etecentexperinrents include 4-wave nrixing spectroscopy of 



single semiconductor quantum dots [12| or nonlinear optical measurements on single elec- 

n n 

trons in F centers[13|]. CARS microscopy [14j is being extended to small interaction volumes. 
The traditional macroscopic formulation does not apply for small samples and a quantum 
electro dynamical (QED) formulation is called for. 

A QED treatment is essential for predicting observables related to the quantum nature 
of the field. A powerful quantum master equation approach has been successfully used in 
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quantum optics to describe fluctuations of the laser field, photon statistics, entanglement 



and squeezing 



151 ]- 25 1 . In this paper we consider simple observables that are usually cal- 



culated by the semiclassical approach. However a quantum description of the fields allows 
a more precise definition of the signal, provides new insights and shows the limitations of 
the semiclassical approach. Our final results are expressed in terms of multipoint correla- 
tion functions of the material system. These can be evaluated for more complicated models 
than the Bloch-type equations ordinarily used in quantum optics. Details of molecular com- 
plexity, coupling to a solvent with arbitrary time scale and excitonic effects may be readily 
described by these expressions jsj]. 

The present approach is based on the many-body non-equilibrium Green's functions 
(NEGF) technique [26I] - 



32j. It provides a single microscopic definition of the signal that 
includes both spontaneous and stimulated emission and generates both incoherent and coher- 
ent signals. We shall use a superoperator in Liouville space representation that is commonly 
applied to study electrical conduction in quantum junctions [331]- [35(. In Sec. [IT] we give 
a general definition of the signal (Eq. (111ft ) and express it in terms of Green's functions 
of the system and the field. The formalism will be demonstrated by first considering the 
signal generated by a single molecule. Applications are made to two types of incoherent 
experiments, spontaneous light emission (SLE) (Sec. IIIII) and to pump-probe spectroscopy 
(Sec. HVj) . A purely diagrammatic derivation of the final expressions for the SLE signal using 
the technique of Keldysh Schwinger loops [3^] is developed. We present rules that directly 
translate the diagrams to partially-time-ordered expressions for the signal, both in the time- 
and in the frequency- domain. This representation, has been shown to be most useful for 
the study of frequency-domain experiments within the semiclassical approximation 36|. We 
extend it to quantum fields, and show that it yields more compact expressions than the 
fully-time-ordered double-sided Feynman diagrams. 

In Sec. |Vl we turn to the response of an assembly of noninteracting particles. The 
expressions of Sec. [Ill are generalized by simply replacing the dipole operator of a single 
molecule by that of an assembly of particles. The signal is separated into an incoherent and 
a coherent part. The latter, which is described by semiclassical optical susceptibilities, is 
dominant in macroscopic samples. However for small iV both have comparable magnitudes 
and coexist. The coherent signal does not exist in the limit of a single molecule. This 
indicates a breakdown of the semiclassical theory which does predict such a signal. Finally, in 
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Sec. I VI I we show that within QED heterodyne detection emerges as an incoherent stimulated 
emission process in the detection mode. This is in contrast with semiclassical treatment 
where we first generate a coherent macroscopic signal, interfere it with a local oscillator and 
then obtain the signal. All of these steps are avoided by the microscopic QED definition 
given here. Sec. IVIII provides a summary of our results. 

II. INCOHERENT OPTICAL SIGNALS 

We first consider a single molecule driven by an optical field E(t). This is all we need for 
calculating incoherent signals. The extension to coherent signals from molecular ensembles 
will be done in Sec. [V] 

The total Hamiltonian is given by 

H = H + H F + H int , (1) 
where Hq represents the free molecule, and 

H f = 2j f^ s a\a s , (2) 

s 

is the radiation field Hamiltonian. In Eq. (j2J) a s (a|) denotes the destruction (creation) oper- 
ator for the sth field mode and a s (a\) satisfy the boson commutation relation [a s , a s ,] = 5 SS >. 
For clarity, we shall omit the unit vector describing the field polarization of a particular 
mode, and treat the field as a scalar. The generalization to include the full tensorial expres- 
sions is straightforward. 

We shall treat the emitted photon modes with frequency u s quantum mechanically. £(£') 
are the positive (negative) frequency components of the quantized electric field 

E{r,t)=£(T,t)+&(r,t) , (3) 

with 

£(r,f) = W^Y\e^^-', (4) 

s ^ ' 

and Q is the quantization volume. 

In the rotating wave approximation (RWA) the molecule-field interaction is given by 

H M (t)=£(r,t)Vi + &(T,t)V , (5) 
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where we have partitioned the dipole operator jx as fi = V + . Here V\V) are the creation 
(annihilation) operators for excitations. The time-dependence in Eq. (JSJ) is in the interaction 
picture with respect to H F . Consequently, the time dependence of the total (field+molecule) 
density operator pr(t) is governed by the Hamiltonian 

H(t) = H + H mt (t) ■ (6) 
We assume that the field is initially in a coherent state, 

|tfjr) = A exp j^a^A |0> , (7) 

where 

A = exp(^|a s | 2 ) , (8) 

s 

is the normalization such that {9f\^f) — 1- in Eq. ([7]) a s is the eigenvalue of a s , a s \^fp) — 
a s \^F), and |0) is the vacuum state. The expectation value of the field is then 

(tt F |£(r,t)|tt F ) = £(r,t) + c.c. , (9) 

where 

S(r,t) = ^(^) 1/2 ^ r -^a s , (10) 

s ^ ' 

is the field amplitude at space-point r. 

We shall define the time and frequency resolved signal S(t) as the rate of change of the 
photon occupation, i.e. 

*<«> == s Wo, T • ^ 

where 

A^:=^aia s , (12) 

s 

denotes the photon number operator. This definition applies to both spontaneous and 
stimulated process, as will be demonstrated later. In Eq. ffl2l the sum extends over all 
modes of interest in the experiment under consideration. 

We shall adopt of the following convention for the various ensemble averages, one of which 
has been used in Eq. (ITTj) . In Eq. (TTTT) . (. . .) D T stands for a trace with respect to the density 



operator of the molecule where Hi nt includes all field modes ("total" system). According 
to the definition of the signal, Eqs. fTTTl) and (fT2j) . the field modes can be partitioned into 
detected (or signal) modes and the incoming modes. For later use, we introduce the average 
(. . .) D , which denotes a trace with respect to the density operator of the molecule calculated 
in the interaction picture where Hi nt only includes the incoming modes. Averages with 
respect to a noninteracting system, where Hi nt =0, will be denoted by (...). This corresponds 
to the limit t — > — oo (adiabatic switching). 

The expectation value in Eq. (TlTI) may be conveniently evaluated by switching to the 
Heisenberg picture, using the basic identity 




and the canonical commutation relations for a s (aj). In Eq. ffl3l) we denote operators in the 
Heisenberg picture by a subscript "H" . 

Evaluating the commutator in Eq. ( Tl3l) and transforming back to the original interaction 
picure, we finally obtain for the signal, Eq. (flTl) . 

5(t) = -^Im{(^(r,t)Vt) DT } . (14) 

Taking into account the time dependence of the total (system + field) density operator, 
which is determined by the Hamiltonian of Eq. (Q, the trace on the rhs of Eq. (Ti^|) can 
be computed perturbatively in the field E(r, t) by switching to the interaction picture with 
respect to the Hamiltonian H . We shall compute this trace using superoperators in Liouville 
space. To this end we associate with each operator A in Hilbert space, a left (L) and a right 
(R) operation in Liouville space, by defining 

A L X := Ax , 

A R X := XA , (15) 

where X denotes an arbitrary Hilbert space operator. Furthermore, we introduce the fol- 
lowing linear combinations of L/R operations, which will be referred to as +/— operations 

A ± := -L[A L ±A B ] . (16) 
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A key tool in the following manipulations is the time-ordering operator in Liouville space, 
T; when acting on a product of superoperators, it reorders them such that their time 
arguments increase from right to left. In the interaction picture Eq. (1141) can be expressed 
as 



T£ L (r,t)v£(t)exp j-~ J drv / 2^-(r)| 



o ( rt 

S(t) = --Im 



(17) 



where 'Hint- is the superoperator corresponding to Hint, Eq. 

Note that H int contains two interaction terms, one with the incoming field modes H' int 
and the other with the signal modes H" nt . As before, the explicit time dependence represents 
the interaction picture 

a. x (t) := ei HFxt a sX e-i HFxt , 

V x (t) ■= e * noxt V x e-K noxt , (18) 

where X = L,R and Hfx (Hqx) is the superoperator corresponding to H F (H ). 
To lowest order in H" nt Eq. (|T7|) gives 

x f dr (ra sL (t)Vl(t)nl t _(T)) (19) 

J -co ' " 

The superoperator expression corresponding to Eq. (jSJ) gives 

V2H int -(T) = £ L (r,T)Vl(r)+£l(r,T)V L (r) 

- vI{t)£ r {t,t)-V r (t)£ r {v,t) . (20) 

Substituting in Eq. ( TT9l) and factorizing the correlation functions into products of field and 
system parts, we get 

S(t) = ^Re^ ss ,(^,) 1/2 e j(ks - k ^- r 
x t)L> m (t, *) - ^(t, r)D LL (r, t)] , 

(21) 

where we have introduced the field and system nonequilibrium Green's functions, dxy and 

Dxy, 

dxY(t,r) := ~(ra sX {t)a\, Y {r))^ , 

D XY {t,r) := -- (rVx{t)V^{r))^ . (22) 



In Eq. ( !2T|) we assume that the system is in its ground state for t — > — oo, hence the following 
system correlation functions vanish 

(^v2(*)Vj(r)) D = 0, 

(rv3(t)Vj(r)) D = 0, f>r. (23) 

From the basic definitions of the Green's functions we have 

D LL {r,t) = e{T-t)D RL {r,t) + e{t-r)D LR {T,t) , 
<kL(*,r) = 9(t-T)d RL (t,T)+9{r-t)d LR (t,r) . 



Substituting Eq. (124"]) in pTl) gives 

S(t) = ^Re^ ss ,(^) 1/2 e i(ks - k ^- r 



x dr [df R (t, t)D rl (t, t) - dg L (t, t)D lr (t, t 



Using Eqs. ([7j) and (122]) . the Green's functions for the field are given by 

dl s R (t,r) = -- (a;,(r)a s (t)) D 



= -~ (a s (t)al,(r)) 



D 



The signal, Eq. (|25|) . then becomes 



= ^i m ^(^ s ,)i/V( k - k ^ r e-^* 



x 



(24) 



(25) 



~e-^*e^' r « s «:, , (26) 



-~e-^V^(<^+a s a:0 • (27) 



/ dre^' r [a,a^D M (r,0-(^ + a a a;)D L «(r,*)l . (28) 

J — oo 



In Eq. ( I2"8"j) . the term with S ss / represents spontaneous emission, whereas the terms pro- 
portional to a s a*i correspond to stimulated processes. Note also that the Green's functions 
Dxy contain the density matrix of the driven system which involves all interactions with 
the incoming field. 
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Equation ([28]) can also be written as 



S(t) = 21m dr [D rl {t, t) - D LR (r, t)] S(t, t)6*(r, r) 
T.s too dre-^-^DLRir, t), 



(29) 



or equivalently 



S{t) = 2lmJ t _ oo dr [D RL (r,t)£(r,t)£*(v,r) - D LR {t,r)£*(r,t)£{r,r)} 



(30) 



In Eq. (1301) we made use of the symmetry relation 



D 



(r,t) 



D LR (t,r) , 



(31) 



to arrive at an expression, where interactions at time t always occur from the left (i.e. on 
the ket). This choice will be adopted in the diagrammatic representation, which will be 
introduced in Sec. IIHI 

Equations (|29l) or equivalently (1301) constitute our basic NEGF expression for incoherent 
optical signals. The first term in the brackets describes the creation of excitations in the 
system by absorbing photons from the field and the second term represents the reverse 
process. Both processes are stimulated. The signal is thus given by the net photon flux. 
The last term represents spontaneous emission. In the coming two sections we shall apply 
these expressions to compute different signals. 

III. DIAGRAMMATIC REPRESENTATION FOR SPONTANEOUS LIGHT 
EMISSION (SLE) 

To describe spontaneous light emission (SLE), which is an incoherent process, we assume 
that the scattered field is initially in its vacuum state and is generated by spontaneous 
emission. Thus £{r,t) = and Eq. ( 130~1) . reduces to 




(32) 
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where it is furthermore assumed that the scattered field has only one mode uj s thus dropping 
the sum over s (signal modes). By integrating over modes we can replace by the density 
of modes. 

Expanding the Green's function D^r to second order in the incoming field yields 



In order to go any further, we need to expand the four-point correlation function of Eq. 
( l33l) in L/R operations. A systematic, diagrammatic technique for selecting the terms in 
Eq. ( 1331) that contribute within the RWA, will be introduced next. 

We start with Eq. ( 1331) . where the correlation function originates from expanding the 
Green's function D^it, r) in Eq. fl32|) to second order in the incoming field. Hence in 
what follows, we associate the time variables T\ and r 2 with the incoming field, whereas 
the variables t and r correspond to the signal field. We also note that thanks to the time 
ordering operator, Eq. (J33l) is symmetric with respect to interchanging the dummy variables 
Ti and r 2 . 

In order to derive an explicit expression for this correlation function in terms of superop- 
erators, we first note that the interactions at times t and r both represent photon emission, 
one on the ket at time t, the other on the bra at time r. Since the system is assumed to be 
initially in its ground state, both emission events have to be preceded by two absorptions 
on either side. Although this implies a certain order of interactions on both the ket and 
the bra individually, their relative ordering in physical time remains a priori undetermined. 
This partial time ordering can be expressed diagrammatically by arranging the interactions 
on a loop. We adopt the following general rules for constructing and reading the diagrams: 

1. Time runs along the loop clockwise from bottom left to bottom right. 

2. The left strand of the loop represents the ket, the right corresponds to the bra. 

3. Each interaction with a field mode is represented by a wavy line on either the right 
(R-operators) or the left (L-operators). 

4. The field is indicated by dressing the wavy lines with arrows, where an arrow pointing 
to the right represents the field annihilation operator S(r,t), which involves the term 



Ssle(£) 



Qh 2 



cu s Re fl^dr f dn / dT 2 e iuJ ^-^ 




(33) 
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e i(k s -r-cj s t) ^ see Q). Conversely, an arrow pointing to the left corresponds to 
the field creation operator £t(r, t), being associated with e -*( k s-r-w)_ This j s made 
explicit by adding the wavevectors ±k s to the arrows. 

5. Within the RWA (Eq. ([5])), each interaction with S(r,t) is accompanied by applying 
the operator V\ which leads to excitation of the state represented by ket and deex- 
citation of the state represented by the bra, respectively. Arrows pointing "inwards" 
(i.e. pointing to the right on the ket and to the left on the bra) consequently cause 
absorption of a photon by exciting the system, whereas arrows pointing "outwards" 
(i.e. pointing to the left on the bra and to the right on the ket) represent deexciting 
the system by photon emission. 

6. The interaction at the observation time t, is fixed and is always the last. As a conven- 
tion, it is chosen to occur from the left. This can always be achieved by a reflection of 
all interactions through the center line between the ket and the bra, which corresponds 
to taking the complex conjugate of the original correlation function. 

7. Interactions within each strand are time-ordered, but interactions on different strands 
are not. Each loop can be further decomposed into several fully-time-ordered diagrams 
(double sided Feynman diagrams). These can be generated from the loop by simply 
shifting the arrows along each strand, thus changing their position relative to the 
interactions on the other strand. Each of these relative positions then gives rise to a 
particular fully time-ordered diagram. 

8. The overall sign of the correlation function is given by (—1) Nr , where Nr stands for 
the number of interactions from the right, at times T\ and t 2 . 

We note that loop diagrams drawn according to the rules presented above, lead to double 
sided Feynman diagrams that follow the standard conventions employed within the semi- 
classical theory of nonlinear optics jsj]. Using these rules, the SLE is represented by the 
single loop diagram displayed in Fig. [U We denote the incoming field as ±ki and the signal 
field as ±k 2 . The loop translates Eq. fl33l) into the following expression for the signal, 

W*) = r^^Re J dr J dr a J dr 2 e^^ 
(TV L (t)^(r)V L t (r 2 )V R (r 1 ))(T^(r,r 2 )4(r,r 1 )) , (34) 
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where an additional prefactor of 2 was added to take into account the symmetry with respect 
to the dummy variables T\ and r 2 . 

If desired, this can further be decomposed into fully time-ordered terms using the double 
sided diagrams of Fig. [2J As indicated above we can now e.g. shift the arrow in Fig. [T] 
corresponding to time r 2 along the ket, thus obtaining 3 possible relative positions with 
respect to the interactions at times r and T\. The loop diagram consequently splits into 3 
double sided Feynman diagrams. These are depicted in Fig. [2j 

The resulting fully-time-ordered expression for the signal reads 

Atx t 

Ssle(*) = Re dr / dn / dr 2 ^*(r, n)£i(r, r 2 ) 

e ^(t-r) [eitrMmMw) {v L {t)vU,r)V R {r 1 )Vl{T 2 ) 
+0(tr 2 )^(r 2 r)^(rr 1 )(v L (t)y2(r2)V^(r)V«(r 1 )^ 

+e(tr)e(TT 2 )e(T 2 T 1 )(v L (t)v^T)vl(T 2 )v R (T 1 ))} . (35) 



Here, we have made use of Eq. ([7j) to evaluate the correlation functions for the incoming 
field, since 



2nh 

Efanfair,^) . (36) 



Note that, in Eq. (}3"5l) full time-ordering is expressed by adding a product of step functions 
to each of the contributing terms. We apply the short hand notation 9{t\t 2 ) = 9{t\ — r 2 ). 

Next, we calculate the signal for a frequency-domain experiment with stationary beams. 
This will also illustrate how to extend the diagrammatic rules to frequency-domain experi- 
ments. We adopt the following convention for the Fourier transform of a function / 



f(cu) = J dte^f(t) , 

f(t) = ^ J da,e--7» • (37) 

When the incoming field is stationary, the integrals in Eq. ( 1341) can be evaluated directly, 
which by application of Eq. ( )3Tj) leads to a multiple Fourier transform of the partially 
ordered 4-point system correlation function in Eq. (134]) . To simplify this calculation, it is 
advantageous to switch to a new set of time variables ("s- variables"), si := t — r 2 , s 2 := t — r, 
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S3 '■— t — Ti, which represent the time intervals of interactions along the loop. This choice 
becomes intuitive when looking at the diagram, Fig. HJ 

We can now show that the correlation function in Eq. (1351) depends on the s-variables 
in a simple way. To this end, we transform the system correlation function in Eq. f[3~4|) to a 
form that only involves L-superoperators, hence 



Mri)Vl(r)V L (t)Vl(r 2 )) (38) 



In terms of the s-variables, Eq. fl38l) becomes 



y L G^{sz)Vfo{s 2 )V L g{ Sl )Vl) , (39) 

where we define the Liouville space propagator 

Q{s) := #(s)exp (--V^Ho-S^ . (40) 

Note that, when associating all interactions with the ket, one needs to propergate twice 
backwards in time. This amounts to applying the retarded propagator twice. 

Using Eq. ( 1391) the integration in Eq. ( 134T) . leads to a multiple Fourier transform of the 
correlation function in terms of the loop (s)- variables. This is readily calculated, giving the 
following expression for the frequency-domain SLE signal 

SsLE^l, UJ 2 ) = 1^1,0 1 w 2 



Im (v L g\u x )Vlg\u x - u 2 )V L G{ui)Vl 



(41) 



Since all interactions in Eq. (T4T1) are "left" -operations that act on the ket, it is possible 
to express the final result for the signal in Hilbert space, giving 

47T 

'S'SLE^I, C0 2 ) = ^r-|£l,o| ^2 

Im (vG\u g + u^G^ujg + u x - u 2 )VG(u g + uji)V^ . 

(42) 
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Upon conversion of Eq. (14T]) . the frequency-domain Liouville space propagators, 

^ ^ hw — \^2Ho- +ir) ' ^ 
are replaced by Hilbert-space propagators, 

6(u> + u> 9 ) = — — * . (44) 

ftcj + — H + t7] 

Here, u; g is the material frequency of the ground state, which accounts for the free evolution 
of the bra. In both Eqs. (1431) and (T44]) the infinitesimal rj > arises from causality and 
guarantees the convergence of the Fourier transform. 
Equation (14*21) can be recast in the form 

a i \ 2ni \c 12 

(45) 

which by using the level scheme of Fig. [3] (a) yields the Kramers-Heisenberg formula 

Q l \ ^ \C 121 ,2 

<JSLEl^l)^2j — —"^7^^21^1,01 |X 



(46) 



where 



Xca = W ; = • (47) 

In Eqs. (146]) and (14~7)) cj ca = aj c — c<j a denotes the transition frequency between the levels c 
and a and we have added a phenomenological dephasing rate r ca . 

By proceeding along the same line for an arbitrary loop diagram, one can establish the 
following rules, that allow to translate a given diagram into its frequency-domain expression. 
These complement the rules given earlier for time-domain expressions. 

1. In the frequency-domain the loop translates into an alternating product of interactions 
(arrows) and periods of free evolutions (vertical solid lines) along the loop. 

2. Since the loop time goes clockwise along the loop, periods of free evolution on the 
left amount to propagating forward in real time (G), whereas evolution on the right 
corresponds to backward propergation (&). 
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3. Each G adds the multiplicative factor 
by {-%). 



i, whereas each & results in a multiplication 



4. Frequency arguments of each propagator are cumulative, i.e. they are given by the 
sum of all "earlier" interactions along the loop. The ground state frequency uj g must 
be added to all Green's functions. 

Equation (l4"2j) can be immediately generated from Fig. @] by applying these rules. 

IV. THE PUMP-PROBE SIGNAL 

The simplest nonlinear optical technique involves two fields, ki (the pump) and k2 (the 
probe). The signal is defined as the difference of the probe transmitted intensity with and 
without the pump. We assume that the probe intensity is high so that spontaneous emission 
is negligible compared to stimulated emission. The second term in Eqs. (1291) and (130]) . which 
describes spontaneous emission, is thus neglected. Using Eq. (|29|) and expanding to second 
order in the incoming field, the signal then reduces to 





(48) 



Using the identity 




(49) 



we can express the pump-probe signal in terms of +/— operators as 




(50) 



This is equivalent to the classical expression for the signal in terms of x ■ This becomes 
even more apparent if Eq. (157J1) is rewritten in the form 




(51) 
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which is equivalent to Eq. (1501) . if only terms proportional to the field intensities are taken 
into account. The details are given in appendix El 

The forgoing example of SLE illustrated the use of Keldysh-Schwinger loops to compute 
the signal. In particular, by using this partially time ordered diagrammatic representation, 
we were able to reduce the number of diagrams from three to one. The merits of this compact 
notation will become even more apparent for the pump-probe signal, as will be shown next. 

We start with an expression for the signal using Eq. (130]) 



where in analogy to Eq. ( )48|) . we neglect the spontaneous emission term. 

We note that the second term in Eq. fl52l) coincides with the SLE correlation function, 
Eq. (I33I) . What remains therefore is to resolve the first term in terms of L/R operations. 
For simplicity we assume a three level ladder with sequential transition dipole moments as 
shown in Fig. [3] (b). 

We again start our analysis by looking at the two interactions at times t and r. Con- 
trary to the second term in Eq. ( 1521) . both of these correspond to absorptions rather than 
emissions. For the dipole selection rules shown in Fig. [3](b), this leaves two possibilities for 
placing the remaining interactions along the loop: We either have three interactions on one 
side, or two interactions on either side (see Fig. H]). 

For the former case and the model of Fig. [3] (b), these can only correspond to two 
absorbtions and one emission. By applying the same argument used for SLE, the earliest 
of these three interactions must necessarily be an absorption. To proceed further, one then 
needs to distinguish whether these interactions occur on the ket or the bra. Since the last 
interaction is fixed to time t, the former case correspondingly allows only one possible loop 
diagram, which is displayed in panel (b) of Fig. HI 

When placing the three integrations on the right strand, one has to allow for all possible 
orderings of their associated times r, t\ and t 2 . At this point we reiterate that when drawing 
a loop diagram, even though the relative time ordering of the ket and the bra interactions is 
unspecified, the temporal order within each strand is fixed. Hence having three interactions 
with the bra, generates 4 loop diagrams, which arise upon performing permutions of the 



S PP {t) = 2^3 Re J^dr J dn / dr 2 
S 2 (r,t)S*(r, r) (TV R (r)V^W^-(ri)H mt -(r 2 )) 




(52) 
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arrows on the right side. Panels (e) - (h) of Fig. H] show these diagrams. Note that in this 
case the loop diagrams are actually fully time ordered. 

We now turn to the diagrams with two left and two right interactions. Since the last 
interaction occurs on the ket at time t, for a system starting off in the ground state, the other 
interaction on the ket also must be an absorption. Finally, we note that when calculating 
the correlation functions Eq. (|52|) . one needs to take a trace in the end, which restricts both 
interactions on the bra to be absorptions as well. Arguing along this line leads to the two 
loop diagrams shown in panels (c) and (d) of Fig. HI The two loops again take into account 
the two possible permutations of placing interactions on the bra. 

In summary, the pump-probe signal can be represented by eight loop diagrams, seven 
stem from the first correlation function in Eq. (|52p . whereas the remaining one is SLE-type. 
The corresponding expression for the signal derived using the time-domain rules is given in 
appendix 03j 

For the sake of completness, we shall compare the QED signal with the semiclassical 
result, where the pump-probe signal is calculated as a third order response in the direction 
+ki — ki +k2. in Note that the diagrams in panels (b) - (h) in Fig. HI correspond to +ki — 
\ti — k 2 . This difference can however easily be resolved by taking the complex conjugate of 
these diagrams, which leaves the expression for the signal in Eq. (1B1I) invariant. Pictorially, 
since this amounts to reflecting the arrows through the center line between the ket and 
bra, we recover the classical combination of wavevectors +ki — ki + k 2 . As was the case 
for SLE, we can generate the corresponding fully-time ordered diagrams from the loops 
displayed in Fig. HI The resulting 16 double-sided Feynman diagrams are summarized in 
Fig. [5j In addition, as has been shown in Sec. II I II for SLE, frequency-domain expressions 
follow naturally from the loop diagrams, since the field permutations are already built in. 
Applying the rules given in Sec. II I II to the diagrams for the pump-probe (Fig. HJ), we can 
immediately write down the frequency- domain signal 
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Sp P (u u cu 2 ) = -^|^, | 2 |^,o| 2 Im {- (V&(u g + Ut)V^G^{u g +u t - u 2 )VG{u g + ui)V*) 
+ (vG\uj 9 + uj 2 )V^G{uo 9 +00!- ujVGfag + u{)V^ 
+ (vG\uo g + ujV&iug + + cu 2 )V^G(u g + CJ!)t/t\ 
+ (v&{u g + uj 2 )VG\uj g + u,+ u 2 )V^G{u g + u x )V^ 
+ (vG\uj g + WiJV'tGt^ +CU2 - U2 )V&(u g + u 2 )V^ 
+ (y&(u g + uo 2 )V^G\uo g + lj 2 - uJV&iug + uo 2 )V^ 
+ (y&(u a + cux)VG\u g + ui + uo 2 )V^G{u g + lu 2 )vA 
+ (V&(u g + uo 2 )VG\uo g + u 2 )Vi&(u g + w 2 )M} . 

(53) 

The eight terms correspond respectively to the eight diagrams in Fig. HI 

V. COHERENT VS. INCOHERENT NONLINEAR OPTICAL PROCESSES 

So far we have focused on incoherent processes. We could thus consider a single molecule 
and simply multiply the signal by N in the end. To describe coherent signals we consider 
V in Eq. (j!4p as the dipole moment of a collection of noninteracting molecules at positions 
v a , i.e. 

V = Y,V a S(r-r a ) , (54) 

a 

where V a now denotes the dipole operator of a single molecule. Subsequent interactions of 
the detecting field can take place with the same (a = 0) or with different (a ^ (3) molecules. 
We shall describe an (n + l)-wave mixing experiment, where n incoming modes generate a 
signal in a new direction. When the signal mode (index s) is initially in its vacuum state, 
the spontaneous emission signal can be described by Eq. (1291) modified to include Eq. (1541) . 
We then get 

5 T = > s Re{l + r^/! 00 dr 

J2^ e -i» s (t-r)^ si r -r a} (t^t)^ t)) J , (55) 
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where we have further added a time-integration to describe a frequency-domain experiment 
[37! ] . Equation (155]) is a general expression for spontaneously generated signals. The a = 
ft and a 7^ ft terms in Eq. (1531) give rise to incoherent (Sj) and coherent (Sc) signals, 
respectively. We thus write 

S T (u s ) = N ■ Sfa,) + N(N - 1) • S c (u.) , (56) 

where the incoherent signal is given by 

SVK) = ^Rejl^dt |* dre-*-^) (t^t)^*)) J . 

(57) 

For uncorrelated particles, we have [V^, Vg] = 0. The a ^ ft sum may thus be factorized 
into a and /?, 



^(r)^(t)) c = (K,L(r)) D (^(t)) ] 



^( T )) D (v2(*)) D • (58) 



This gives rise to the coherent term in Eq. (1561) . 

Sofa) = 77^ l^(^)| 2 



»+oo 

P(u s ) = I dTe*Wv2(r)) • (59) 



OO 



D 



For this we have also made use of 



1 



+00 pt 

dt dr = ^ / dt I dr , (60) 



in Eq. 0571) . 

In Eq. (I59p we assumed that the sample is much smaller than the optical wavelength 
or that we have exact phase matching Ak = 0, where Ak is the difference between the 
wavevector of the signal mode and the sum of wavevectors of the n incoming fields. More 
generally, we should replace the N(N — 1) factor in Eq. (1561) by 

F(Ak) = ^2 e~ iAk '( r ™~ r/3 ) . (61) 

a,/3 



For macroscopic samples, the number of terms with a 7^ ft in the double sum of Eq. ( I6TT) 
will by far exceed the terms with equal indices. Evaluating this double sum in the contin- 
uous limit and letting the sample volume go to infinity, gives the standard phase matching 
condition F(Ak) — ► <5(Ak), which yields a directed signal. 
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For small samples both coherent and incoherent terms need to be considered. It is in- 
teresting to note that for a single molecule we never get Sc, we only have Si. Nonlinear 
susceptibilities, even though they are calculated for single molecules, do not represent a spon- 
taneous wave-mixing experiment performed on a single molecule. Using the semiclassical 
approach one may conclude that coherent signals are possible even from a single molecule. 
The QED approach shows that this is not possible. 

Below we compare the expressions for Sj and obtained by expanding Eqs. floTj) and 
(|59|) for n = 1,2, 3, where n denotes the number of incoming modes. We demonstrate that 
the incoherent signal for n incoming modes, relates to a 2(n + l)-point correlation function 
in the dipole moment, compared to an [n + l)-point correlation function contained in p( n \ 

For n — 1 

P (1) K) = fdte^f^dr (TV L (t)Vl(r))e(i,T) , 

gives Rayleigh scattering, whereas 

Sf\u s ) = jdt dre**^) / dn J dr 2 (TV^V^V^n)^)) S*(r, n )S(r, r 2 ) 

is responsible for SLE (Raman and Fluorescence) as discussed in Sec. IIHI 
For n = 2 

P( 2 V S ) = Jdte^ t j t _ oo drjdT l (TV L (t)Vl(r)V-(Ti)\ S(t, r)£*(r, n) 



4l|-[43j, and 



describes 3-wave mixing e.g. second harmonic generation 

Sf (u s ) = Jdt ft^ dre^C-) / dn / dr 2 / dr 3 / dr 4 £ *(r, n)S(r, r 2 )S*(r, r 3 )S(r, n ) 



may represent many possible processes, depending on the level scheme and off-resonant 



detunings, e.g. two photon fluorescence 
For n = 3 

P^(cu s ) = Jdte^J^drJdnJdT^T^ir^S^n) 

represents 4-wave mixing e.g. third harmonic generation and 

Sfius) = jdt dre^^ j dn j dr 2 / dr 3 / dr 4 / dr 5 / dr 6 £*(r, Tl )£(r, r 2 )£*(r, 
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can represent many possible processes, e.g. three photon fluorescence. 

A nanocrystal is made out of many unit cells. Its dipole operator is given by jl — Y2 a A"? 
where fi a represents an electron-hole pair on unit a. It thus behaves as a collection of many 
molecules and it can show spontaneous coherent response 
we have a single particle. 



9] . N is thus large even though 



VI. HETERODYNE-DETECTION AS A STIMULATED WAVE MIXING 

When detecting nonlinear optical signals in macroscopic samples, one can either measure 
the intensity (homodyne detection) or the magnitude and phase of the signal field (hetero- 
dyne detection). The later technique uses an intense external field £lo, usually referred to 
as "local oscillator" , which interferes with the created nonlinear signal in the same direction, 
but is assumed not to interact with the molecule. Using the standard semiclassical approach 
to nonlinear optics, the quantity detected in a heterodyne measurement is |3J] 

/+oo 
dtlm{S* LO (t)P s (k s ,t)} , (62) 
-oo 

where P s is the polarization in the signal direction k s . With the QED formalism developed 
here, we can greatly simplify the derivation of Eq. (1621) and show that it emerges naturally 
as an incoherent rather than a coherent signal. 

Equation (1551) describes the signal generated from an assembly of particles by spontaneous 
emission. Applying Eq. (|54l) to Eq. (|29|) . we can give a similar expression for the stimulated 
signal 

S s Ut) = -fRe/^dr [(rV a , R (r)Vl L (t)) d - (TV a , L (r)V^(t)) ^] S(r,t)S*(v,r) . 

(63) 

In analogy to Sec. (jV]), we can now split Eq. (|63|) into a coherent (a 7^ (3) and an incoherent 
(a = (3) contribution. Since the dipole operators of different, uncorrelated particles a 7^ /3 
commute, the coherent contribution vanishes identically. The present formalism hence yields 
the general result that stimulated processes are always incoherent in nature. In particular 
this applies to heterodyne experiments where the local oscillator serves as the stimulating 
field. It is hence sufficient to look only at a single particle and use Eq. (TH|) derived in Sec. 
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Equation f|T4l) is reminiscent of the expression for heterodyne detection (Eq. (|62l) ). if one 
only looks at one particular mode "s" in the definition of the signal (see Eqs. (fTTj) and (fl2l) ). 
thus 



2 

Sswm(*) = ~^ Im 



T 



(64) 



Here, by adding a subscript a, we make explicit that Eq. (1641) is the signal generated by 
one molecule at position r a . For an assembly of noninteracting particles, this has then to 
be summed over r Q . The only difference of Eq. ( 164"1) as compared to Eq. (1621) . is that it 
contains a field operator S, rather than a classical field £. 

Heterodyne detection can be viewed as an incoherent stimulated emission process in 
the detected mode. Since Eq. ( 1641) is already first order in the mode s, in a subsequent 
perturbative expansion of the density operator the detecting mode will only contribute to 
higher order. Hence, to first order in the coupling with the detecting mode, the time- 
dependence of the density operator is only given by all other modes s' ^ s. The field 
operator therfore acts directly on the state of the system (molecule + field) at t — > — oo, 
which when assuming a coherent state of the field in the same limit (see Eq. ([7])) yields 

S(i) = ^Im[£;(r a ,t)(v)J . (65) 

In Eq. (1651) we additionally took the complex conjugate within the imaginary part. 

The expectation value in Eq. fl65|) can be further expanded in the incoming modes. For 
n modes this gives the classical nth order polarization 

P (n) (r,t) = p(")(t) e <^v s k s')- r . (66) 

Consequently, integrating over all times, we obtain the following generic result for an (n+ 1)- 
wave mixing process (i.e. 1 detected mode "s" + n incoming modes u s' ^ s") 

5(Ak) = - J dflm ft pW(f)/(Ak)] , 

/(Ak) = £V iAk - r * , (67) 

a 

where we explicitly display the spatial dependence of the fields, e.g. 

S s (r,t)=£ s>0 (t)^- r . (68) 
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In Eq. (1671) we have added the sum over r a to make the transition from the signal generated 
by single molecule to one of several molecules (incoherent process). 

Under exact phase-matching conditions (Ak = 0), or if the sample is much smaller 
than the optical wavelength, we have /(Ak) = N. Rather than creating a coherent signal, 
propagating it and interfering as done in the semiclassical approach, we can simply describe 
it microscopically as an incoherent signal. 

Equation fIBTj) . hence allows us to interpret heterodyne detection as follows: ^//fields (the 
n fields (uv, k s >) as well as the detected field (u 8 , k s )) interact with the molecule. This leads 
to stimulated emission of photons of frequency u s . A sum over all molecular positions has to 
be performed (Eq. (1671) ). Going to the continuum limit, which is justified for a macroscopic 
sample, this sum can be replaced by an integral over the sample volume. Extending this 
integral over the entire real space, we obtain the phase- matching condition /(Ak) — > 5(Ak). 
Therefore, in a macroscopic sample the signal is finite only when the stimulating field S s (r, t) 
fulfills this condition. 

VII. CONCLUSIONS 

A microscopic QED treatment of nonlinear optical processes induced by a weak quantized 
field is developed in this paper. In the traditional, widely-used approaches to nonlinear 
optical response, the matter-field interaction is treated semiclassically, hence the signal field 
is obtained by simultaneously solving macroscopic (Maxwell-) and microscopic (quantum 
Liouville-) equations. Here we treat the entire process as a single event. 

Based on a microscopic definition of the spontaneous and stimulated signal, general ex- 
pressions are derived for the nonlinear optical signal involving non-equilibrium Green's func- 
tions in the incoming field. Depending on the experiment under consideration, these can 
then be expanded order by order in the incoming field modes. Application is made to spon- 
taneous light emission and pump-probe spectroscopy. Coherent and incoherent processes 
involving both spontaneous and stimulated emission are treated using a unified framework. 

A diagrammatic derivation of the contributing terms within the RWA using the tech- 
nique of Keldysh Schwinger loop is developed. This respresentation reflects the partially 
time-ordered nature of the Green's functions and hence yields more compact expressions 
than the well established fully time ordered double sided Feynman diagrams. The loop dia- 
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grams are particularly useful for frequency-domain measurements where the bookeeping of 
time ordering is not necessary anymore. If needed for time-domain experiments, the fully 
time ordered expressions can be obtained from the loop diagrams using a simple prescrip- 
tion. Rules for constructing and reading these diagrams are laid out, making it possible to 
intuitively derive the expressions for the signals. 

The practical merits of this partially time-ordered approach are illustrated for the pump- 
probe technique where the number of diagrams is reduced from 16 to 8. This example also 
demonstrates an interesting fundamental merit of the QED treatment of nonlinear optical 
processes, even though in pump-probe spectroscopy quantum effects of the field (i.e. terms 
stemming from spontaneous emission) are usually negligible: Semiclassically, the pump- 
probe signal is obtained by calculating a third order response induced by two fields in the 
direction ki — ki + k 2 . The resulting signal field is then obtained using heterodyne detection, 
where mode 2 acts as its own local oscillator field. Even though the final expression for the 
signal is identical to a QED result, the semiclassical approach creats an artificial asymmetry 
between the local oscillator (mode 2), which by definition does not interact with the molecule, 
and the remaining fields (mode 1 and 2), which give rise to the nonlinear polarization. This 
asymmetry, which is eliminated here, is a direct consequence of the limitations imposed by 
separately treating matter and field within the semiclassical approach. 

Finally, using the QED formalism, we are able to treat both coherent and incoherent 
processes in a unified way. This distinction comes about by calculating the signal for an 
assembly of molecules rather than for a single particle. The signal then splits into two parts, 
a coherent term scaling ~ iV(iV — 1) and an incoherent term ~ N, where N denotes the 
number of particles. As the number of particles in the sample is increased, the coherent 
term becomes more pronounced. For experiments carried out with only a few molecules 
however, both coherent and incoherent terms will generally make comparable contributions 
to the signal. In single molecule experiments only the incoherent term survives. 

Since incoherent signals may not be calculated semiclassically, this leads to a striking 
consequence. For an (n + l)-wave mixing processes carried out on a single molecule we show 
that the signal is not given by an (n + l)-point response function of the dipole operator 
related to P^ n \ as predicted by the semiclassical approach. It is rather given by a different 
2(n + l)-point combination of correlation functions. This doubling arises since the signal 
may not be recast as an amplitude square and we must calculate the signal itself, not an 
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amplitude. The present formalism allows a microscopic calculation of nonlinear optical 
experiments on single molecules which have only become feasible recently. 

Finally, we address this apparent limitation of the semiclassical description from a more 
general viewpoint and show that heterodyne detection for an (n + l)-wave mixing experiment 
can be simply viewed as an incoherent stimulated emission process in the detected mode. In 
contrast, homodyne-detected n-wave mixing is a coherent spontaneous emission process. For 
a macroscopic sample, by summing over the molecules in the interaction volume, we recover 
the phase-matching condition, i.e. the heterodyne-signal is finite only if the stimulating 
mode (i.e. the local oscillator) matches this condition. 



APPENDIX A: DERIVATION OF EQ. (j5lj) 

In this appendix we illustrate the equivalence of the QED expression for the pump-probe 
signal with the classical signal expressed in terms of x • Starting from Eq. ( 1501) and 
applying the definition of Eq. ( fT6l) . we can write 

1 



+ 4(r,T)y_ + £i(r,T)yJ . (Al) 



Making use of Eq. flAll) . the correlation function in Eq. fl50l) can be recast in terms of 
correlations of the system and the field. Calculating the field correlations, only 6 of the 
possible 16 terms are nonzero, giving 



2 (TK 1 r (t)\/_(r)H' mt _(r 1 )^' mt _(r 2 

+ ^TV5(t)y_(r)VL(ri)yi(r 2 )) (r^(r,n)4(r,T 2 

+ (rv5(t)VL(r)Vi(n)y + (r 2 )) (rf+(r,ri)fi(r,r 2 

+ (TVl(t)V4T)V + ( Tl )vl(T 2 )) (T£Ur,r 1 )£ + (r,T2 

+ (TVlWV^V-intflfo)) (r^(r,r a )^-(r,r 2 

+ (TVl(t)V4T)Vl(n)V„(T 2 )) (T£4r, n )£l(r,T 2 ) 



(A2) 
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where the system is assumed to be in its ground state for t — > — oo, hence 

(TVlW^TW-^V-^)) = . (A3) 

Note that since Eq. (J5D|) is symmetric with respect to the variables T\ and r 2 , Eq. flA2[) 
simplifies further and only includes three distinct terms upon integration, i.e. 

/ dn / dr 2 2(T\/ R (r)V r i t (t)H / m4 _(r 1 )H' mt _(r 2 )) = 
2 / dn / dr 2 x 

+ (r^(i)y_(T)y_(r 1 )^(T 2 ))(r4(r,T 1 )^_(r,T 2 )) . 

(A4) 

For a classical incoming field, the correlation function in Eq. (1A2j) is equal to only the first 
term on the rhs, where 

(T£ + (r, Tl )£l(r,T 2 )) = e(r, Tl )S*(r,T 2 ) + ^e^"^ . (A5) 

For a classical field the second term in Eq. (1A5I) will be negligible. 

The last two terms in Eq. (1A4I) are due to the quantum character of the incoming field, 
since 

(r^+Cr.n^Kr.ra)) = ^e^^^n - r 2 ) , 

(r4(r,rO£-(r,T 2 )) = -^e^^^n - r 2 ) . (A6) 

Hence, neglecting the last two terms in Eq. ()A4j) and making use of Eq. (|A5j) . we obtain 
Eq. (ED). 

APPENDIX B: TIME DOMAIN EXPRESSIONS FOR THE PUMP-PROBE SIG- 
NAL 

In this appendix we give the expressions for the pump-probe signal in the time-domain 
obtained by applying the rules summarized in Sec. IHII to the loop diagrams in Fig. HJ This 
complements the frequency-domain expression given in Eq. (|53|) . Applying the time-domain 
rules for the loop diagrams in Fig. HJ we obtain 
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S PP (t) = F Re {iL ^ I / dr 2 5 2 * (r, t)5 2 (r, r) 
Ty x (t)Vj(r)^ t (r 2 )y i? (r 1 ))(r4 i (r,r 2 )4 t iii (r,r 1 ) 
+ j! 00 dr/dri/dr 2 ^(r,t)£J(r,r) 



+ 



TVl(t)V R (T)V L (n)Vl(T 2 )) {TSUv, rO^Cr, r 2 ) 
Yy x t (t)y Ji (r)yi(r 1 )y ii (r 2 )}(r^ i (r ) r 1 )4 t iii (r,r 2 ) 
+ (rvJWVkWyJC^ViiCra)) (r^fon)^ (r,r a ) 



(Bl) 

where we kept the order of the diagrams in Fig. |H Similar to Eq. ([Mil , an additional factor 
2 accounts for the symmetry with respect to T\ and r 2 . This can again be broken into fully 
time-ordered terms using the double sided Feynman diagrams displayed in Fig. [5j We then 
get 

Spp(t) = ^Re y* dr y dn y dT 2 £ 2 *(r,i)£ 2 (r,-r) 
£*(r, t 1 )S 1 (t,t 2 ) [9(^)9(^)9^) {Mt)V ] R (r)V R (ri)Vl(r 2 )) 
+ 0(f7s)0(T 2 T)0(T7i) (^(O^^)^^)^^)) 

+ 0(tr 2 )0(r 2 r)0(m) (^(t)^^)^^)^^)) + 0(tr)0(r^)0(r 2 ri) (^(t)^)^^)^^) 

+ £ 2 (r,t)£ 2 *(r,r) {5 1 *(r,n)5 1 (r,r 2 ) [0(tT)0(<rri)0(7i7*O (^(^^(r)^^)^^)) 

+ 0(tri)0(rir)0(rr 2 ) (^(^^(tOVr^)^^)) + 0(tri)0(rir 2 )0(r 2 r) (^(^^(ri)^^)^^ 

+ 0(tri)0(rir)0(rr 2 ) (^(^^(tOFrCt)^)) + 0(tr)0(rr 2 )0(r 2 ri) (^(^(t^t,)^) 
+ ^^^(nr) (v2(*)Vj(r 2 )y2(n)V^(r)) +0(tr 1 )0(r 1 r 2 )0(r 2 r) (^(^(n)^)^ 
+ 0(tr 2 )0(r 2 r)0(m) (^(^(^(t)^) 



+ £ 1 (r,Ti)£J , (r,T 2 ) [0(t7-)0(rri)0(rir 2 ) (^(^(t^ti)^ 

+ 9(tr 2 )9(r 2 r 1 )9(T 1 r) (vl(t)M'ni)Vkn)V R ('r)) + e(tn)9(TiT)e(TT 2 ) (vl(t)Vl(Ti)V R (r)V R (T 2 



+ 9{tn)9(r 1 r 2 )9(r 2 r) (y L t (t)V^(r 1 )V R (r 2 )V R ( 



(B2) 



Comparison of Eqs. (1B1I) and (1B2[) illustrates the tremendous simplification achieved by 
the more compact partially-time-ordered approach. Employing the diagrammatic Keldysh- 
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Schwinger loop technique, the number of diagrams is reduced by half (16 vs. 8 diagrams). 
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FIG. 1: Loop diagram for SLE. ki is the incoming field, and k 2 is the signal field. Note that the 
interactions are time ordered within each strand, but not between strands. s\, s 2 , S3 are the delay 
times between the interactions along the loop. 
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FIG. 2: Time-ordered Feynman diagrams of SLE, generated by shifting the arrows of Fig. Q] along 
each strand thus changing their relative time ordering. Each of the possible three relative positions 
then gives one fully time ordered diagram (double sided Feynman diagram). 



31 



/'/» 



M«6 



(a) (b) 

FIG. 3: The three level systems with sequential dipole couplings used in the derivation of the 
Kramers Heisenberg relation for the SLE (panel (a); Eq. (|47p ) and the pump-probe signal (panel 
(b))- 
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FIG. 4: The eight loop diagrams for the pump-probe signal (Eq. ([53]) ) . Note that diagram (a) 
coincides with the SLE (Fig. [I]). 
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FIG. 5: The Double sided Feynman diagrams resulting from Fig. [Hupon breaking the loops into 
time ordered contributions. Panel (a) - (d) are the time ordered diagrams corresponding to the 
loops shown in panel (a) - (d) of Fig. [H respectively. Since the four loops with three interactions 
on the bra (Fig. HI panels (e) - (h)) are already fully-time ordered, each gives a single double sided 
diagram. These are not repeated here. Overall, the eight loop diagrams of Fig. Oyield 16 Feynman 
diagrams. 
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